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Abstract 

A three-state model based on the Potts model is proposed to simulate financial mar- 
kets. The three states are assigned to "buy", "sell" and "inactive" states. The model 
shows the main stylized facts observed in the financial market: fat-tailed distributions 
of returns and long time correlations in the absolute returns. At low inactivity rate, the 
model effectively reduces to the two-state model of Bornholdt and shows similar results to 
the Bornholdt model. As the inactivity increases, we observe the exponential distributions 
of returns. 

Keywords: Econophysics; financial market; Potts model; fat-tailed distribution 

1 Introduction 

Recent empirical studies of financial prices such as stock indexes, exchange rates estab- 
lished the stylized facts for price returns [Q: fat tailed distributions, absence of autocorre- 
lations in return, long time correlation in absolute return, etc. A variety of models have 
been proposed to imitate the real financial market and to understand the origins of the 
stylized facts|H]-|12j. Many models are able to capture some of the stylized facts. Among 
them, simple models may serve to find the universality behind the financial market 9 . 

Bornholdt proposed a particularly simple spin model for financial markets based on 
the Ising model [T2]. His model has been further investigated by several authors [131 1141 
I15[ I16|. In the Bornholdt model the two states of the Ising spins are assigned to "buy" or 
"sell" orders. The essential ingredients of the model are two conflicting interactions: (i) 
nearest neighbor interaction (ferromagnetic) (ii) global interaction (anti- ferromagnetic) . 
The nearest neighbor interaction gives an effect to imitate neighbors. On the other hand, 
the global interaction coupled to the magnetization plays a role of the minority effect. 
These interactions cause a complicated dynamics, i.e., there is no single stable phase, 
rather the model exhibits metastable phases. As a result, the distribution of the price 
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returns shows a fat-tailed distribution and the volatility or the absolute return is strongly 
clustered which results in the long time correlation. 

Motivated by the two-state Bornholdt model, we propose a three-state model based 
on the Potts model. Compared to the Bornholdt model, the model includes an additional 
state which may be assigned to "inactive" state. Therefore the model incorporates three 
states assigned to "buy" , "sell" and "inactive" . We simulate the model with the additional 
inactive state and study whether the results show the stylized facts observed in the real 
financial market. 



2 Potts-like Model 

Following the heat bath dynamics proposed by Bornholdt jT2] . we construct a model based 
on the three-state Potts model. Each agent i lives on a site of (i-dimensional lattice and 
has a spin Si. First let us consider the Bornholdt model. Spins Si have two states (+1, — 1) 
and are updated to 5- with the following probabilities 1 . 

P{Sl -> S^ = eMP(h itS , - aSiSi\M\))/Z (1) 

where 

Ks>. = S S jt S' t > (2) 

<ij> 

and the summation is taken over the nearest neighbors j of the site i. Z is the 

normalization factor determined such that the following equation is satisfied. 

£ P( Sl ^k) = l. (3) 
fc=i,-i 

M is the magnetization given by 

1 N 

M = ]V^>' (4) 

i 

where N is the total number of agents. 

The factors Ph ijS > represent the nearest neighbor interactions according to which agents 
imitate their neighbor. On the other hand, the factors aSiS'^M] give the minority effects. 
Namely if the magnetization \M\ gets large, the agents tend to change their states. 

Next we add an additional state corresponding to the inactive one. We assign "0" to 
this state. We define the probability to have this inactive state as 

P(S l - Si = 0) = exp{(3(h iiS , - e ia \M\ - 2 7 K))/Z, (5) 

where K is the inactivity rate given by 

1 N 



K = ^J2 S s„o, (6) 



1 Here we use the two-state Potts representation rather than the Ising one. 
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and ej = 25s u o ~ 1 which takes 1 for Si = and —1 for Si ^ 0. The factor ej«|M| 
corresponds to the minority effects and 2^fK controls the inactivity rate. By varying 7 we 
study the model for various inactivity rates. 

In summary, our updating scheme for the three-state model is as follows. We update 
spins Si to S'i with the following probabilities given as 

P( a J expMhs'-aSiSUMDyz for S> = 1,-1 

1 ij \ exp((3(h iiS < -6ia\M\- 2yK))/Z for S' { = lU 

where Z is the normalization factor determined by 

£ P(5i^fe) = l. (8) 

fc=l,-l,0 



3 Results 

In our simulations, we use a 100 x 100 square lattice with the periodic boundary condition. 
When we update spins, we choose a site randomly and update it with the probabilities 
given by eq.Q. After 100 2 chosen sites are updated, time is incremented and we also call 
this one sweep. The simulations are started with a random configuration. The first 10 4 
sweeps are discarded as thermalization and we accumulate 10 6 sweeps for analysis. We 
use the random number generator in Ref.|17j for a single processor of SGI Altix3700. 

The simulations were done at j3 = 4 and a = 10. This value of j3 is larger than the 
critical coupling (3 C ~ 1.005[T5] of the three-state Potts model on a 2-dimensional lattice. 
In other words, the temperature T = 1/(3 is below the critical temperature. We define the 

1 N 

return through the mag netization[T5] as r(t) = M(t+1)-M(t) where M(t) = ^ X) 5 *^)' 

i 

In our model 7 controls the activity rate of agents. Fig.l shows the average inactivity 
rate (K) as a function of 7. As 7 increases, the inactivity decreases. For 7 > 1 the 
inactivity rate is very small: for 7 > 1 the inactivity rate < 10 -3 . Thus for 7 > 1 the 
model is expected to effectively reduce to the two-state model similar to the Bornholdt 
model. 

Fig. 2 shows the distributions of returns at j3 = 4 and a = 10 for various 7. The 
distributions are not the Gaussian distributions but shows fat tailed ones. The shape of 
the distribution at 7 = 1 is similar to that observed previously in the Bornholdt mo del [To]. 
As 7 decreases, however, the distribution changes its shape. At 7 ~ 0.05 the shape seems 
to be an exponential. We will study the shape of the distribution in detail with the 
cumulative distributions. 

Fig. 3 shows the cumulative distributions of the absolute returns \r\. For 7 = 1.0, 
the power-law behavior ~ r~ v is seen in a range of 0.002 < r < 0.1. The exponent v 
is estimated to be a = 2.3 which is the same value obtained in Ref.^SJ. This is not 
surprising because at this 7 the model is effectively a two-state model. 

As 7 decreases, the power-law behavior disappears and the shape of the distribu- 
tion seems to change to the exponential at large r. Assuming the exponential form 
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Figure 1: Average inactivity rate at j3 = 4 and a = 10 as a function of 7. 



exp(— 7]r/a r ), we obtain 7/ = 1.07 at 7 = 0.05 for r > 0.004. a r is the standard deviation 
of r, a 2 = ((r — (r)) 2 }. Recent studies claim that the distributions of returns in emerging 
markets show the exponential distribution. The values of r\ are evaluated to be 77 ~ 1.3 
and 1.5 in the Indian market49 a , and 1.7 in the Brazilian market 20 . Furthermore there 
is also a claim that at mesoscopic time scale the distributions in developed markets also 
show the exponential distributions 21 . 

The autocorrelation function C(t) of a time series g(k) is defined as 



C(t) 



{g{k)g(k+t))-{g{k)f 



of, 



(9) 



where a 2 = (g(k) 2 ) — (g(k)) 2 . Fig. 4 shows the autocorrelation functions C(t) of returns. 
As observed in the real financial market, the autocorrelation of returns disappears quickly 
as the time increases. An interesting observation in Fig. 4 is that the autocorrelation is 
negative at small t. This negative autocorrelation is also seen in the financial market fQ. 

Contrary to the autocorrelation of the returns, the autocorrelation of the absolute 
returns shows long time correlation. Fig. 5 shows the autocorrelation of the absolute 
returns for 7 = 1.0 and 0.05. At 7 = 1.0, we do not see a power-law behavior but 
rather observe an exponential behavior. Assuming exp(— t/to), we obtain to ~ 1280 or 
l/t w 7.8 x 10~ 4 by a fit in a range of 400 < t < 3000. Thus the autocorrelation 
function is a very slowly decaying exponential. The similar exponential behavior of the 
autocorrelation function is also seen in the Bornholdt model [Ej. As 7 decreases, the 
autocorrelation starts to show a power-law behavior. At 7 = 0.05 we observe C{t) ~ t~ 
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Figure 2: Histogram of returns at (3 — 4 and a = 10 for various 7. Each histogram except at 
7 = 1.0 is shifted up adequately to avoid overlapping each other. 
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Figure 3: Cumulative histograms for 7 = 1.0,0.5,0.1,0.08,0.05 from bottom to up. A line 
proportional to r -2 3 and a curve of exp(— 1.07r/a r ) are also drawn. 
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Figure 4: Autocorrelation function of the returns for 7 = 1.0 and 0.05. 




4 Summary 



We proposed a three-state model based on the Potts model. The results shows the main 
stylized features: fat-tailed distributions of returns, long time correlations in the absolute 
returns. 

At low inactivity rate (7 > 1) the model effectively reduces to the two-state model 
similar to the Bornholdt model. We observed a power-law behavior for the cumulative 
distribution of the returns: ~ r -2 ' 3 . The autocorrelation of the absolute returns does not 
show a power-law behavior but an exponential behavior. 

As the inactivity increases, the power-law behavior for the cumulative distribution of 
the returns disappears and the cumulative distribution shows an exponential behavior. 
Although the empirical study of the financial markets shows a power-law 22 a ~ r -3 , the 
exponential behavior observed in our model is not peculiar. Recent studies [TU1 I20j claim 
that emerging markets such as those in Brazil and India show the exponential distributions 
of returns. Moreover, the distributions of the developed markets at mesoscopic time scale 
also show exponential behavior 21 . 
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